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Abstract. In this study we explore, both from theoretical and experi- 
mental side, the effect of Fe doping in ZrC>2 (ZrO-z'-Fe). By means of 
first principles simulation we study the magnetization density and the 
magnetic interaction between Fe atoms. We also consider how this is 
affected by the presence of oxygen vacancies and compare our findings 
with models based on impurity band |T] and carrier mediated magnetic 
interaction [2] . Experimentally thin films (~ 20nm) of ZrO^ '-Fe at high 
doping concentration are grown by atomic layer deposition (ALD). We 
provide experimental evidence that Fe is uniformly distributed in the 
ZrOi by transmission electron microscopy (TEM) and energy disper- 
sive X-ray (EDX) mapping, while X-ray diffraction (XRD) evidences 
the presence of the fluorite crystal structure. Alternating gradient force 
magnetometer (AGFM) measurments show magnetic signal at room 
temperature, however with low magnetic moment per atom. Results 
from experimental measures and theoretical simulations are compared. 

1 Introduction 

Dilute magnetic semiconductors (DMS) are materials in which magnetic impurities 
are introduced in order to produce a magnetic ground state. These systems have 
received great attention in recent years, since the discovery of carrier induced ferro- 
magnetism in (In, Mn)As [3] and (Ga, Mn)As [4], and are believed to be fundamental 
to fabricate spin-based electronic devices. Recently a new class of DMS has been in- 
vestigated, namely DMS based on high-fc oxides, i.e. dilute magnetic oxides (DMO), 
after the experimental evidence of room temperature magnetism in transition met- 
als (TMs) doped zirconia |5l6l7l8j (Zr0 2 ), hafnia [giTU] (Hf0 2 ), and titania [TTTT2j 
(Ti0 2 ) and the theoretical prediction of high T c in TMs doped Zr0 2 |13I14] . 

The understanding of DMS /DMO physical properties constitutes a challenge for 
the theory as the fundamental mechanism leading to ferromagnetic (FM) interaction 
between the dopants cannot be explained in terms of simple exchange mechanisms, at 
least at low doping concentration, being the latter often too short-ranged. Among the 
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other, two theretical models have been proposed to describe FM effects: the first is 
based on the presence of impurity states in the crystal, impurity band model (IBM) [1] , 
the other on carriers in spin polarized bands, carrier mediated model (CMM) [5] which 
is indeed a refined version of the Zener model. 

From the experimental side the inclusion and influence of magnetic dopant, such 
as Fe, Co, Ni and Mn, is not clearly understood. Indeed, while several DMS/DMO 
have been predicted to have a Curie temperature (T c ) above room temperature, no 
experimental report of T c > 300K has been left unchallenged by other studies [IS] . 
Moreover some results suggest that magnetic impurities act as paramagnetic (PM) 
centers with unusually long relaxation time |16) . at least at very low doping concen- 
trations. 

In this manuscript we study iron doped zirconia (ZrC>2'-Fe) focusing our atten- 
tion on the magnetic properties of the system. In Sec. [2] we provide a structural 
characterization of thin films grown by ALD. In particular we show that the dop- 
ing is uniform, with high iron concentration and no segregation, and that zirconia is 
in the tetragonal/cubic structure. Thus theoretically we focus our attention on the 
tetragonal structure of zirconia with substitutional iron doping uniformly distributed 
in the sample. In Sec. [3] we study the magnetization of the system and how it is 
influenced by defects, i.e. oxygen vacancies (Vb), comparing our results with the IBM 
and the CMM. Indeed in a recent work we showed that iron doping induces Vb, with 
a ratio yv /Fe — 0.5, for charge compensation, and that ZrOi'-Fe films growth by 
ALD presents a ratio close to one half [T7]. We finally investigate the magnetization 
of the films growth by ALD. From the magnetization at saturation we extract the 
magnetic moment per atom which is discussed in view of the results from theoretical 
simulations. 



2 Structural Properties of Zr0 2 :Fe 

In order to describe the effect of iron doping in zirconia, ZrO%:Fe thin films were 
grown on Si/Si02 substrates in a flow-type hot wall atomic layer deposition reactor 
(ASM F120) starting from /3-diketonates metalorganic precursors. Ozone was used as 
oxidizing gas in the reaction process. The Fe concentration in ZrO^'-Fe films can be 
tuned tailoring the Zr/Fe precursors pulsing ratio. In the present work however we 
focused our attention on the high doping regime keeping the pulsing ratio fixed. The 
growth temperature was maintained at 350°C. After the deposition the films were 
annealed at 600°C in N2 flux for 60s to study their thermal stability. Further details 
on the samples preparation can be found in Ref. |18j . 

Film crystallinity was checked by X-ray diffraction (XRD) at fixed grazing in- 
cidence angle lu = 1° and using Cu K a (A = 0.154 nm) monochromated and colli- 
mated X-ray beam (Italstructure XRD 3000, details on the measurements can be 
found in Ref. |19]V In the present work all measures shown are from the same 
film which we chose as representative of the high-doping concentration samples. 
Fig [TJa) shows that the films present a cubic/tetragonal crystalline structure with 
an estimated cell parameter, assuming a cubic cell, a = 5.024 A with a contrac- 
tion of about (a exp — a® xp )/a® xp w —1.0%; here a? exp = 5.074 is the lattice pa- 
rameter measured for undoped ZrOi films. Theoretically we found, at x — 25%, 
{dtheo — a theo) / a theo ~ ~ 0.4%, with a® heo — 5.11, the lattice parameter computed for 
Zr0 2 , while the tetragonal deformation reduces from 3.05% to < 1.%. 

The samples were then characterized by TEM, performed by a JEOL 2200FS 
microscope equipped with a high-angle annular-dark-field (HAADF) detector, in- 
column energy filter and EDX spectrometer. The layered geometry of the samples 
is shown in fig. [T{6). The Fe atomic concentration measured by EDX is 25%, in 
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Fig. 1. (colors online) Strucutral properties of iron doped zirconia. Panel (a): XRD patterns 
show that the film presents the cubic/tetragonal phase |28j . A representation of the fully 
relaxed theoretical structure is also pictured with Fe atoms in red, Zr atoms in blue and 
oxygens in black (we used the Xcrysden software Panel (b): HAADF image and HREM 

inset showing the layered geometry of the sample. The EDX chemical maps show a uniform 
doping distribution with no segregation. 



agreement within few percent to the estimation done by X-ray photo-emission (XPS). 
The EDX maps, taken in different regions of the samples, showed that Fe is uniformly 
distributed across the film and no clusters or segregation at the grain boundaries have 
been found (see Fig.QJfr), insets). 

Thus theoretically we computed, from first-principles, the ground state of the 
the tetragonal phase of ZrC>2'-Fe with uniform iron doping distribution. We used 
the PWSCF [50] package considering a super-cell with 96 atoms; for all systems the 
atomic positions are fully relaxed [21] . The ground state was computed with the gen- 
eralized gradient approximation [22] (GGA) to the density functional theory (DFT) 
scheme 23 24J with ultra-soft pseudo-potentials 25 26 . We have recently shown that 
iron, in the zirconia host lattice, behave like Yttrium (V) which is among the most 
studied dopant of this oxide: it replaces zirconium atoms in the ZrC>2 lattice inducing 
Vo for charge compensation, with a ratio yv /Fe = 0.5, and stabilizes the tetrago- 
nal phase above xf £ ~ 11% (27117] , Calculations have been done at different doping 
concentrations, XFe = 6.24, 12.5, 25.0% and different Vo concentration, considering 
the Vo to Fe ratios yv /Fe — 0, 0.5, 1, focusing our attention at the high doping 
concentration limit, i.e. the concentration measured experimentally. In the inset of 
Fig. a relaxed structure at yv /Fe = 0.5 and xp e — 25%, which is the expected 
structure in our films [17], is shown. 



3 Magnetic Properties of Zr0 2 :Fe 

Starting from the fully relaxed structures considered at XFe — 25% we computed 
the local magnetization of the system. In Fig. [2] the magnetization density, i.e. the 
difference between the spin-majority and the spin-minority density is represented 
for yvo/Fe = 0, 0.5, 1. In order to compare the present results with the IBM also 
the total density of states (DOS) are plotted. The FM configuration is considered for 
better clarity as it helps to distinguish between the minority and the majority spin 
channel. However all the conclusion we will draw in the following also hold for the 
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Fig. 2. (colors online) Theoretical density of states (panels a-c) and magnetization density 
isosurfaces at m = 0.01 a.u. (panels d-f) for the ferromagnetic configuration at yv /Fe = 
0, 0.5, 1. The density of states projected on the d atomic orbitals is also shown. The crystal 
field splitting is visible with the e doublet, d z i and d x 2_ y 2 in blue, at higher (lower) energy 
in the majority (minority) spin channel respect to the t triplet, d xy , d xz and d yz in red. In 
panels d-f Fe atoms are represented in red, Zr atoms in blue and oxygens in black. 



anti-FM configuration (see Fig.[3J&-c) for a comparison of the anti-FM magnetization 
density). 

At yvo/Fe = 0, the system shows holes in the valence band (Fig. E^a)) because 
iron act as an acceptor; the projected DOS shows that these are localized on the 
Fe(d xy ) level. As a consequence the majority-spin d band is not completely filled, 
the magnetic moment per iron atom is 4 \ib and iron is forced in a +4 oxidation 
state. This also breaks the spherical symmetry of the system and the shape of the 
magnetization density close to the oxygens is not isotropic, with a peculiar doughnut 
shape oriented in the xy-plane. 

Creating Vo, i.e. increasing yv /Fe > 0, electrons are released in the system. 
These fill the empty d xy levels with a charge-transfer mechanism. At yv /Fe = 0.5 
the majority d band is completely filled, the magnetic moment per iron atom is 5 hb 
and the magnetization has a spherical shape. The system turns into a charge-transfer 
semi-conductor (Fig. Hffe)). At yv /Fe = 1 electrons fill the minority d levels and the 
magnetic moment per iron atom is reduced to 4 fig. However the minority d electron 
does not participate in the bonds, thus the anisotropy is only weakly transferred to 
the p orbitals of oxygen (Fig.^c)). In all configurations the magnetization is mainly 
located around the Fe atoms and it is in part transferred to the oxygen next nearest 
neighbor suggesting a short range magnetic interaction (Fig. W[d-f)). 

It is worth to compare the DOS of our system with the IBM. A crucial assumption 
of the IBM is that the presence of Vo creates impurity states with poorly localized 
electrons. These, already at low doping concentration, would overlap to create an 



Will be inserted by the editor 



5 



impurity band that can mediate the FM interaction among nearby iron atoms. Indeed 
Vo are a common defect in pure ZrOi films where they create impurity states close 
to the conduction band, however in the case for ZrOi'-Fe the situation is different, as 
shown in Fig. ^a-c). For low Vo concentration, i.e. yv /Fe < 0.5 no impurity states 
are associated to Vo while for x > 0.5 impurity states appear but higher in energy 
than the empty spin minority (i-levels. Thus the extra electrons are trapped in the 
Fe(d) levels and the impurity states remain empty. More in general, the majority 
of magnetic transition metals (TMs), i.e. Fe, Co, Ni, Mn, and Cr, have +2 or +3 
as most stable oxidation state and we expect a similar picture for all TMs doped 
XO2 oxides, at least for yv /Fe < 0.5 (1.0), in the case of +3 (+2) oxidation state. 
Thus, according to our results, the IBM cannot be invoked to describe the magnetic 
ground-state for Zr02-Fe and for X02-TAI in general. 

The situation depicted is much closer to the CMM proposed for Ga\- x Mn x As. In 
Ga\- x Mn x As Mn acts as an acceptor which compensates the antisite defects com- 
monly present in GaAs giving a charge-transfer semiconductor; in ZrO^'-Fe Fe acts 
as an acceptor compensating the Vo ■ The main difference between the two systems is 
the nature of the host lattice, a covalent semiconductor the former and a polar oxide 
the latter. In Ga\- x Mn x As, when XMn exceeds the antisite defect concentration, the 
system behaves like a metal with holes in the valence band. In Zr02 : Fe instead, at 
UVo/Fe < 0.5, the holes locate onto the d xy states (see Fig.^fr)); conduction could be 
possibly obtained only at high doping concentration with a hopping like mechanism. 

We finally focus our attention on more standard exchange mechanisms studying 
the energy of the magnetic interaction in our system. This in principle can be com- 
puted as the energy difference between the FM and the PM phase. The latter however 
can be hardly described within a periodic code, as the description of random magnetic 
moment orientations would require huge supercells with a non-collinear description 
of the wave-functions. The energy difference per iron atom between the FM and the 
anti-FM configuration however can be used as a reasonable approximation [15] . We 
found that at lower doping, xp e = 6.25, 12.5%, the energy difference is very low with 
AE/kp of the order of few kelvin [30]. The anti-FM configuration is slightly favored 
at yvo/Fe — 0.5, while the FM one is slightly favored at yv /Fe — 0. Instead at 
XFe = 25% in both cases the anti-FM configuration is favored with AE/kp ~ 15QK 
for yvo/Fe — 0.5 and AE/kp ~ 20K for yv /Fe — 0. This suggests a short range in- 
teraction which becomes relevant at xf s > xp, with xp the percolation threshold, i.e. 
an anti-FM super-exchange mechanism which is dominant in the yv /Fe = 0.5 case. 
For the yv /Fe = configuration however this super-exchange mechanism appears 
to be in competition with a FM interaction which could be explained in terms of the 
CMM model, though with a weaker effect due to the low mobility of the holes. This 
results in a weakly FM interaction at low doping and a weakly anti-FM interaction 
at higher doping. 

The magnetic properties of sample with xp e « 25% were studied by means of 
alternated gradient force magnetometry (AGFM). Magnetization measurements, per- 
formed at room temperature by applying the magnetic field parallel to the film plane, 
show a clear hysteresis loop characterized by a coercive field uqH c « 0.03 T (Fig. 
G2a)). From the saturation magnetization value a magnetic moment per Fe atom 
m « Q.Qup can be extrapolated. This is consistent with the values reported in the 
literature but is one tenth of the theoretically predicted value, m ~ 5^ip. Moreover 
theoretically we found that the anti-FM configuration is the most stable at least for 
yvo/Fe = 0.5 which is the expected situation in our films |T7|. In Fig.^b-c) we plot- 
ted the magnetization density in the anti-FM configuration for yv /Fe = 0.5. As in 
the case of FM ground state, the magnetization is mainly located on the Fe atoms 
and the next nearest neighbor oxygens. 
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Fig. 3. (colors online) Panel a: magnetization hysteresis loop measured by AGFM at room 
temperature. The resulting magnetization at saturation can be interpreted in terms of a 
anti-ferromagnetic structure if yv Q /Fe deviates from 0.5 of about 10%. The magnetization 
density for the anti-FM configuration at yv /Fe = 0.5 is plotted in panels (b)-(c). Panel b: 
2D contour plot of the magnetization density along the plane (1 -1 0). Panel c: Positive (red) 
and negative (blue) isosurfaces of the magnetization density at m = ±0.02. 



The discrepancy between the theoretically and the experimental results could be 
possibly explained supposing that experimentally yv /Fe = 0.5 were not fulfilled by 
few percent of the Fe atoms, which thus would be in the Fe A+ /Fe 2+ oxidation state. 
If we assume that these ions would couple with FM interaction, with a CMM like 
mechanism, the net results would be a weak magnetic system with low magnetic 
moment per atom. However further investigations are needed, both experimentally 
and theoretically, to clarify this point. 



4 Conclusions 



In conclusion we have studied both theoretically and experimentally thin films of iron 
doped ZrC>2- Experimentally we have shown that iron distributes uniformly at high 
doping concentration with no cluster formation and that the system is magnetized at 
room temperature. With first principles simulations we have discussed the possible 
magnetization mechanisms comparing our results with other model proposed in the 
literature. We showed that the impurity band model cannot be invoked in the case 
of iron doped zirconia while the carrier mediated model could be possibly considered 
for uncompensated oxygen vacancies concentration (i.e. yv /Fe 0-5). The standard 
super-exchange mechanism however appears to be dominant for yv /Fe = 0.5. 
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From magnetization measurements of highly doped samples we found a low satu- 
ration magnetization value corresponding to low magnetic moment per Fe atom. We 
tried to interpret the experimental results in view of our theoretical findings. However 
the difference between theory and experiments suggest that other effects need to be 
taken into account for a correct description of the magnetic properties of the system. 
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